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Abstract 

I attempt to construct U(l) conformal field theory(CFT) in the Tomonaga-Luttinger(TL) 
liquid with 1 jr^ long-range interaction (LRI) . Treating the long-range forward scattering as 
a perturbation and applying CFT to it, I derive the finite size scalings which depend on the 
power of the LRI. The obtained finite size scalings give the nontrivial behaviours when (3 is 
odd and is close to 2. I find the consistency between the analytical arguments and numerical 
results in the finite size scaling of energy. 

1 Introduction 

Electron systems have attracted our much attention in the low energy physics. As the dimension 
of the electron systems decrease, the charge screening effects become less important. In spite of 
these facts, models with short-range interaction have been adopted in many researches of one 
dimensional electron systems. The recent advanced technology makes it possible to fabricate 
quasi-one-dimension systems. Actually in low temperature the effect of Coulomb force has been 
observed in GaAs quantum wires [Tj , quasi-one-dimensional conductors O El H] and ID Carbon 
nanotubes 0I51E]. 

The systems with 1/r Coulomb repulsive forward scattering was investigated on the long 
distance properties by bosonization techniques [SJ. The charge correlation function decays with 
the distance as exp(— const. (In x) 1 / 2 ) more slowly than any power law. The momentum distri- 
bution function and the density of state does not show the simple power law singular behaviour. 
The logarithmic behaviours appear in the power [9]. These mean that the system is driven to 
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the Wigner crystal which is quite different from the ordinary TL liquid. The investigation for 
the interaction l/r 1 ~ e through the path integral approach [10] reconfirms the slower decaying of 
the single particle Green function for e = and leads the faster decay for < e(<C 1) than any 
power type. 

The numerical calculation in the electron system with the Coulomb interaction shows that 
the larger range of the interaction causes the insulator (charge density wave) to metal (metallic 
Wigner crystal) transition [11]. In the spinless fermion system, the convergence of the Luttinger 
parameters exhibits the quasi-metallic behaviour different from the simple TL one |12j . 

As I will discuss below, the forward scattering is irrelevant for (3 > 1. As an instance of the 
effect of the long-range Umklapp scattering, it was reported that the 1/r 2 interaction makes the 
system gapless to gapful through the generalized Kosterlitz-Thouless transition [13j . 

In this paper I discuss CFT in the system with LRI. The basic assumptions of CFT are 
symmetries of translation, rotation, scale and special conformal transformation. Besides them I 
assume the short-range interaction in the CFT. Hence it is a subtle problem whether the CFT 
can describe the system with LRI. 

Of LRIs, up to now, the solvable models with 1/r 2 interaction were discussed |14 4 115 1 [TB I 117]. 
With the Bethe ansatz, the conformal anomaly and the conformal dimensions were calculated 
and the system proved to be described by c = 1 CFT. In fact the central charge from the specific 
heat agrees with c = 1. On the other hand, the ground state energy is affected by the LRI and 
the periodic nature. The effective central charged deviates from c = 1. 

In general, the CFT for LRI which breaks the locality, has been left as unsettled problem. It 
is significant to clarify the validity of the CFT to the systems with LRI. I investigate the tight- 
binding model with l/r@ interaction as one of such problems. The low energy effective model 
consists of TL liquid, the long-range forward scattering and the long-range spatially oscillating 
Umklapp scattering. Extending arguments appearing in Ref. [18] to the TL liquid with the 
long-range forward scattering, I derive the finite size scalings. In the tight-binding model with 
\j T P interaction, I calculate numerically the size dependences of energy and the coefficients of 
1/L y . And I see numerically the relations between the velocity, susceptibility and Drude weight, 
which CFT requires. 

2 Field theoretical approach 

I consider the following tight-binding Hamiltonian of the interacting spinless Fermions: 

L L 

h = - j>W + m + 1 2> - uma - *> - 1/2), w 

X I define the "effective central charge" by c = — in the finite size scaling of the ground state energy E g — 
aL — . I use the word "effective central charge" in this sense. 
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where the operator cj (cJ) annihilates (creates) the spinless Fermion in the site j and pj = CjCj 
is the density operator. In order to treat this model under the periodic boundary condition, 
I define the chord distance between the sites i and j: r%j = (— sin nfcj) ) where L is the site 
number. Using this, I express the LRI as V(i — j) = — T ^u-iW „ • 

(_ sin " L JU )P 

By the bosonization technique, I obtain the effective action of the Hamiltonian (pQ) for the 
arbitrary filling: 

S = drdx — (V(j)) 2 +g / drdxdx d x (f>(x , r)V (x — x )di<j)(x ,r) 

J 27TA J 

+ g j drdxdx cos(2fcpx + V2<f>(x, r))V(x — x ) cos(2fcpx + V2<ft(x , t)) , (2) 

where V(x) = K is the TL parameter and kF is the Fermi wave number. And g is the 
couping constant proportional to g. The first term of ([2]) is the TL liquid and the second term 
is the long-range forward scattering. The last term is the spatially oscillating Umklapp process 
which includes cos 2y2(f) which comes from the interaction between the neighbour sites. 

Schulz analyzed the effects of the Coulomb forward scattering by the bosonization technique 
in the electron system|8j. He discussed the quasi- Wigner crystal of electrons due to the Coulomb 
forward scattering. Here I focus on the effects of the 1/r 13 forward scattering in the spinless 
Fermions system. I treat the action : 

drdx- — — (V0) 2 + g j drdxdx d x (f>(x, t)V(\x — x \)d x /<p(x , r) (3) 



2-kK 

for any filling hp. To investigate in the Fourier space, I choose the form V(x) = ^ 2+ ^ 2 ^/ 2 , 
where a is the ultra-violet cut-off. In the wave number space, the action ([3]) is expressed as 

/2ir 
dqdw{-(q 2 + w 2 ) + gq 2 V Xq)}\<Kq, w)\ 2 , (4) 

where V{q) is the Fourier transformation of V(x): 

n«) = r( ^ /2 2 ) ^ F / 2 -i /2 (^) /3/2 ' 1/2 ^-v^ a g)- ( 5 ) 

Here K u {x) is the modified Bessel function of z^th order and T(x) is the gamma function. From 
this, the dispersion relation is 

W " = 9 2 {1 + ^^)}- (6) 
The long wavelength behaviors of V(q) are given by 

A + B(qa) 2 + C{qaf- 1 + ■■■ (3 > and f3 + odd 

A + B In ga H = 1 

V(q) ~ { A + B{qa) 2 \nqa + C{qa) 2 + ■■■ (3 = 3 (7) 

A + B(qa) 2 + C{qaf ln(ga) + D{qa) A + ■■■ (3 = 5 
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where A, B, C and D are the functions of j3. For the case where (3 > and ^ odd, the coefficient 
B = B((3), C = C{p) is given by 

vr 3 / 2 f 
B{J3 ^ ~ 4 2 / 3 / 2 - 1 /2r(^)r(/3/2)sin^^ 

c <« - -' 3/2 ^- lr( m )r ; /2)sin i^ ' < 8 > 

From the eqs. ([6]) and ([7]), I see that (</a)^ _1 and In go terms for < (3 < 1 affect the linear 
dispersion essentially. Especially for (3 = 1, there is the analysis by Schulz, where the charge 
density correlation function is calculated [8]. According to it, in the present spinless case, the 
LRI drives the ground state from the TL liquid to the quasi- Wigner crystal as (3 — > 1+. The 
slowest decaying part of the density correlation function is given by 



< p(x)p(0) > ~ cos(2/cpx)exp(— cy^og x), (9) 

where c is a function of K, which exhibits slower spatial decay than the power decay of TL 
liquid. 

Then I see the effects of the long-range forward scattering in the standpoints of the renor- 
malization of g. The renormalization group eqs. of g, v and K are simply derived for the 
long wave-length (see Appendixes.). From the renormalization eqs., the g terms are relevant for 
(3 < 1, marginal for (3 = 1 and irrelevant for (3 > 1. Thus it is expected that the system becomes 
the quasi- Wigner crystal caused by the forward scattering for (3 < 1 and the system becomes 
the TL liquid when (3 > 1. I see that 

= J dx'd x <j)(x,T)V(x -x')d x ><t>(x ,t) (10) 

has the scaling dimension x g = (3 + 1 for 1 < (3 < 3 and 4 for (3 > 3. As the weak logarithmic 
corrections appear for /3=odd, I here distinguish <3?(x) for /3=3 from the scaling functions. I also 
find the consistency on these scaling dimensions by CFT. By using the first order perturbation, I 
can know the effects of the long-range forward scattering. Based on CFT, the finite size scalings 
of energies for no perturbations are given [191 1201 EH [22] by 

„ ^ 2irvx n 

= — 

Eg = e 9 L -jr> (11) 
where x n is the scaling dimension of the primary field denoted by n, v is the sound velocity and 
c is the central charge. Considering the LRI, I can extract the corrections to these energy size 
scalings(see Appendixes.): 

AE n = ^(l + ^^ + 0(l/L 2 )) 
E g = ( eg+9 (0)const.)L-— (c + 5 (0) const. +5(0)-^^ + 0(1/L 2 )), 

(12) 
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where /?(> 1) is not odd. And the constants are the functions of (3. Note that for (3 = odd 
cases, the logarithmic corrections appear. They correspond to the integer points of the modified 
Bessel function, which appear in the long- wave behaviours ©• I can reproduce these anomalies 
for (3 = odd by the CFT. Moreover from CFT I can show that there are the anomalies in the 
general excitations and the ground state energy. The details are shown in Appendixes. The 
0(1/L 2 ) terms come from the irrelevant field L_2-^-2l and the long-range g term. The first eq. 
of (|12p means that the long-range forward scattering has the scaling dimension x g = (3 + 1 
for 1 < (3 < 3 and 4 for (3 > 3 effectively. These respective scaling dimensions are consistent 
with the estimation from the renormalization group eqs. of g, that I mentioned above (see 
Appendixes.). 

The energy finite size scalings (|12p mean that the LRI has the higher order influences than 
1/L to the excitation energy and the LRI affects the 1/L term in the finite size scaling of the 
ground state energy. Here I note that it becomes difficult to calculate the central charge from 
finite size scalings ([lip unless the effects of the LRI to 0(1/ L) terms are known. 

It is notable to compare the eqs. ()12p with the case where the perturbations are of short- 
range type. Ludwig and Cardy calculated the contributions of the short-range perturbation |18j . 
The results for the irrelevant perturbation, — g^2 r 4>(r), which has the scaling dimensions x > 2 
are 



A En = ^(i + ^cw^r 2 ) 

j-j Xfi j-j 

7T11 

E g = (e g + g^) CO nst.)L-—(c + g(Qf^= + 0(l/L^)), (13) 



7™/ /^2 COI1S t- /rlT-fi\ 



where the O(g) terms do not appear in the ground state scaling because we set (<p) = for 
the short-range interaction. These results mean that the x > 2 irrelevant field has influences of 
the higher order to the finite size scalings (jlip . And their result contains parts not so simple. 
There are the special points of scaling dimension x = 1,3,5, and x = 2 which is related to the 
appearance of logarithmic corrections. 

To the contrary, I see ($) ^ in the long-range case, where $ is defined in eq. (jlOp . The 
LRI gives the 0(\/L) intrinsic influence to the finite size scaling of the ground state energy, as 
appearing in the scalings (|12p . even if the LRI is irrelevant, that is, x g > 2. 

3 Numerical calculations 

Through the Jordan- Wigner transformation, I transform the model ([1]) to S = 1/2 spin Hamil- 
tonian for the numerical calculations: 

H = -E( s / 5 7+i +h - c ) + fE 5 ^-^i) 5 i- ( 14 ) 

3 i+3 
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I impose the periodic boundary condition S^+i = Si to this model. Using the Lanczos algorithm 
I perform the numerical calculations for the Hamiltonian (|14p . 

I have found analytically the corrections to the energy scalings (|lip caused by the long- 
range forward scattering. If the oscillating Umklapp process term of ([2]) is irrelevant and does 
not disturb the energy scalings, the finite size corrections due to the forward scattering are 
expected to appear in the excited state energies and the ground state energy. I attempt to 
detect the contribution of the forward scattering. 

I numerically calculate the size dependences of the excitation energy AE(m = 1/L) and the 
ground state energy E g (m = 0), E g (m = 1/L) for g = 0.5. Here I define the magnetization 
m = Y2j Sj/L w hich is the conserved quantity. Fitting the one particle excitation energy as 
LAE(m = 1/L) = a+-£s + jjz,l show the power c versus the powers j3 in Fig. [TJ I see the power 
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Figure 1: The numerically calculated powers c in the excitation energy LAE(m = 1/L) are 
shown versus (3 for g = 0.5. Here I use the scaling form: LAE = a + + where a, b, c and 
d are determined numerically. If the LRI is not present, the energy finite size scaling must take 
the form :LAE = A + -^ + 0(1/L 4 ), where A and B are constant values. 



c agrees with theoretical predictions: (3 — 1 except for (5 = 2. 1 shall discuss the (3 = 2 case later. 
Fitting the ground state energy per site as E g /L = a + jp + , I plot the powers d versus (3 
in Fig. [2j I see that the power d do not show agreements with theoretical predictions (3 + 1 in 
E g (m = 0)/L. These disagreements may be caused by the oscillating Umklapp process which 
becomes relevant at only m = filling. On the contrary, the oscillating Umklapp process is 
irrelevant at m / 0. Actually, in Fig. [51 I see that the power d show agreements with theoretical 
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Figure 2: The numerically calculated powers d in the ground state energies E g {m = 1/L)/L and 
E g (m = 0) / L are shown versus (3 for 51 = 0.5. Here I use the scaling form: E g /L = a + + 
where o, 6, c and d are determined numerically. If the LRI is not present, the energy finite size 
scaling must take the form :E g /L = A + + j^, where A, B and C are constant values. 
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predictions (3+1 in E g (m = 1/L)/L except for (5 = 2. 

As I stated above, for (3 = 2, the power c in the excitation energy LAE(m = 1/L) = 
a + jc + -rz apparently shows disagreement with theoretical value (3—1 and likewise for (3 = 2, 
the power d in the ground state energy E g (m = 1/L)/L = a + + jd apparently shows 
disagreement with theoretical value (3 + 1. I investigate the reason for these disagreements. 
In Fig. [3] I show the numerically obtained coefficient of 1/L d in the size scalings E g (ra = 
0)/L, E g {m = 1/L)/L = a+ + and the numerically obtained coefficient of 1/L C in the size 
scaling LAE(m = l/L)=a + j^ + jj I . I observe that the coefficient of 1/L C in LAE(m = 1/L) 
and the coefficient of 1/L d in E g (m = 1/L)/L become small around (3 = 2. So for (3 = 2, 
1/L 2 dependence appears rather than 1/L in LAEim = 1/L) (see Fig. [TJ). Likely for (3 = 2, 
1/L 4 dependence appears rather than 1/L 3 in E g (m = 1/L)/L (see Fig. [2j). I observe that the 
coefficient of 1/L d in E g (m = 0)/L show the different behaviour from those in E g (m = 1/L)/L 
in Fig. This difference may come from the spatially oscillating Umklapp process that opens 
the gap at m = and disturbs the finite size scaling. 
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Figure 3: I show the numerically obtained coefficients of 1/L d in the size scalings E g (m = 
0)/L, E g (m = 1/L)/L = a + + j-& and the numerically obtained coefficient of 1/L C in the size 
scaling LAE(m = l/L) = a+^ + -^. I observe that the coefficient of 1/L d in E g (m = 1/L)/L 
and the coefficient of 1/L C in LAE(m = 1/L) become small around (3 = 2. The coefficients of 
1/L d in E g {m = 0)/L show the different behaviour from that in E g {m = 1/ L)/L. This difference 
may be caused by the spatially oscillating Umklapp process term. 



I can obtain A{(3) in the scalings (|40p and (|45p by evaluating the integrals. The results are 
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shown in Fig. 0] (a) and (b). The analytical A{(5) in the scalings (|40p and A(f3, s) in the scalings 
(|45p for s = fit with the points in Fig. [3J well. The curve only for s = in Fig. IHb) shows 
the good fitting. This point shall be discussed later. These reveal that the present numerical 
calculation of the tight-binding model agrees with the CFT analysis of the long-range forward 
scattering. 




3.5 




Figure 4: (a) A(P), the coefficient of in the eq. (I40p is shown. I see that A(j3) has zero 

point close to j3 = 2. This curve coincides with the results from the numerical calculation in the 
tight-binding model shown in Fig. [3j (b) A((3), the coefficient of 1/L@, in the eq. (|45p is shown 
for some s. Analytically only s = is meaningful for particle excitations. A(f3) for s = has 
zero point close to (3 = 2. This coincides with the results from the numerical calculation in the 
tight-binding model shown in Fig. [31 

Next I survey whether the long-range tight-binding model satisfies the necessary condition of 
CFT. The operator cos V2<ft has the scaling dimensions K/2 and the operator e =tl v / 2 e h as 
in the regime of the TL liquid. The two quantities 2K/v and vK/2 are the compressibility and 
the Drude weight respectively in the regime of the TL liquid. If c = 1 CFT is valid to the 
tight-binding model with the LRI, the two quantities are related to the two excitations with the 
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symmetries q = ir, m = and q = ir,m = 1/L respectively: 

2K/v = l/(LAE(m=l/L,q = ir))=x 

vK/2 = LAE(q = ir)=D. (15) 

I show the numerically calculated quantities \ and D in Fig. [5] and [U where I use the sizes L = 
16, 18 and 20 and extrapolate the data. For g < 0, x (which is the susceptibility, irrespective 
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Figure 5: The extrapolated K/v(= x/2) is plotted versus the strength g. I use the scaling form 
v/K{L) = v/K(oo) + jjhr for < 3, and v/K{L) = v/K{oo) + ^ for p > 3, where v/K(oo), 
a is determined numerically. 

of the CFT arguments) exhibits the rapid increase which suggests the phase separation. In spin 
variables' language for (pQ), this phase separation is nothing but the ferromagnetic phase. Hence 
for the larger (3 the point of the phase separation approaches to —1. For g > I see the weak 
tendency that the the quantity x becomes smaller as (3 is smaller for g less than about 1. I find 
that the the quantity D of g > become larger as (3 approaches to (3 = 1. 

In Fig. [7]I plot the velocity versus the strength g for the various powers (3, where the velocity 
is defined by 



^-AE(q = 27T/L). 



(16) 



I see that the velocities are finite values for f3 > 1 , as is expected. There are the points where 
the velocities are zero, implying the phase separation. 

In Fig. [8] I plot the quantity versus the strength g for the various powers f3. If the 
present system is described by c = 1 CFT, this quantity is 1 from eqs. (|15p . I find the regions 
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Figure 6: The extrapolated vK(= 2D) is plotted versus the strength g. I use the scaling form 
LAE = a + -jj, where a,b and c are determined numerically. 




g 



Figure 7: The extrapolated spin wave velocity v is plotted versus the strength g. I use the 
scaling form LAE = a + where a,b and c are determined numerically. 
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where = 1 in Fig. [gj The regions become wider as (3 approaches to 1 for g > 0. For larger 
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Figure 8: The normalization is plotted versus the strength g. 



g, the normalization breaks owing to the generations of mass. 



4 Discussion 



I have investigated the system with the 1/r 13 interaction by applying CFT to it and by the 
numerical calculation. At first I have analyzed TL liquid with the l/r@ forward scattering by 
utilizing the CFT and I have found that the l/r@ forward scattering works as higher order 
corrections in the excitation energy, whereas the effective central charge in the scaling of the 
ground state energy depends on the interaction and it deviates from c = 1. The deviation are 
like the solvable 1/r 2 models |14| [15} [TBI fT7] . Next I have numerically calculated the ground 
state energy and excitations energies in the tight-binding model with l/r@ interaction, which is 
expected to include the above \/r@ forward scattering in the low energy. The numerical results 
are in accordance with the analysis with CFT of the long-range forward scattering. Furthermore 
I have numerically checked the normalization -^r = 1, which is the necessary condition for c = 1 



CFT. 



For (3 ~ 2, the coefficient A{f3) in the ground state energy vanishes. This seems to correspond 
to the exact solution for (3 = 2 [17] which states that the finite size scaling of ground state has 
no higher order term than 1/L. The coefficient D((3) of 1/L 3 in eq. (I42p does not vanish for 
(3 = 2. However the present argument is the first order perturbation theory. With higher order 
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treatments, I may clarify this. In any case, with consistency in many points I could construct 
CFT in the system with non-local interaction. 

The numerical calculations in the tight-binding model support the finite size scalings (|40|) 
and (|45p . In one particle excitation energy LAE(m = l/L), the coefficients of l/L^ fit with 
s = case in Fig. SJ The coefficients from the long-range forward scattering are related with 
the operator product expansion. I can prove that only s = case is relevant for the particle 
excitation. Using < ip(z)(p(z') >= — ln(z — z) and < (p(z)(p{z) >= — -f- ln(z — z'), I confirm 
the operator product expansions: 

dtp{z):<*&^: = 1 , : e ^(*W) . +reg . 

4 z — z 

3<p(z) : ^ '>*') : = ^JL. : e*^<-V) . +reg . (17) 



4 


z — z 


iv^ 


1 


4 f 


=- '. 

1 

— z 
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1 


4K (z 


- z'y 


1 


1 


4K (z 


- z'y 



T(z) : e *^(,V) : = : e^V) : : e ^*(*V) :) : +reg. 

iv z — z 

T(z) : e ^(-'- ; ') : = -L 1 : e ^(,',;') . . e ^(,', 2 "') :) . +reg _ ; 

1 '' ' ' K z — z 

where I define T(z) = -j?(d<p(z)) 2 , f{z) = --|(^(z)) 2 and = j?((p(z) - <p{z)). From 

the first and the second eqs., I see C a \o = —iV2/4, C al Q = iy/2/4 for a = 1 and C a \Q = C al Q = 
otherwise, where 0( ) and 1 denote d<p(z) ( d(p{z)) ) and : e l ^ s ( z ' z ) : . From the third and 
the fourth eqs., I see : e l ^ e ( z ^ z ) • have the conformal dimension (1/4K, 1/4K) and spin 0. As 
i(dip(z) — d(p{z))/2 is associated with d a {a) for z = exp(^^), I obtain 



< a\d a (a)\l > 



^i(C 110 -C m )/2 = ^ fora = l 
otherwise, 



which means that only s = is relevant for the particle excitation and the last eq. in (|43p has 
no cosine term. 

I would discuss the size effects for (3 = 1. As seen in the eqs. (|29|) and (|45|) . the velocity shows 
the weak divergence for the size and the Luttinger parameter vanishes gradually for increasing 
size. This is consistent with the numerical tendency (see Figs. 8 and 9 in Ref. [E]). The size 
effect of the Drude weight (proportional to the charge stiffness) is now given by vk/2 ~ const, 
as the logarithmic contributions cancel. The numerical data (see Fig. 7 in Ref. [12]) shows 
the metallic behaviour at small and intermediate magnitude interaction strength (larger than 
the CDW transition point V = 2 by the short range interaction). I think that the long-range 
forward scattering enhances the metallic character. For fairy large interaction the long-range 
Umklapp scattering becomes relevant and the charge stiffness is suppressed. Finally I would 
add the size effect of the compressibility: 

LAE(p = l/2 + l/L) = 1/x = OQnL) -> oo, (19) 
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which comes from the results (|29[) by the RG analysis and the CFT arguments. The compress- 
ibility x goes to weakly for increasing size. 

To summarize, within the perturbation theory I have constructed CFT in the TL liquid 
with l/r@ long-range forward scattering. I have found that the interaction gives the nontrivial 
behaviour for (3 =odd and j3 ~ 2. I have numerically checked the finite size scalings obtained 
from CFT in the tight-binding model with l/r@ LRI. Our analysis and numerical calculations 
exhibit consistency with each other. 

Appendixes 

1. Renormalization group equation 

0</3<lorl</3<3 

I derive the renormalization group equations heuristically. Let us start from the action (J3|): 

S = E E §(9 2 +™ 2 M^)| 2 +5£ E q 2 V(q)\<t>(q,w)\ 2 
w q=—A w q=—A 

A/6 -A/6 A A/6 -A/6 A 

= Ei E + E + Ei+^E{ E + E + E)- ( 2 °) 

w q=-A/b <j=-A g=A/6 w q =-A/b q=-A q=A/b 

The partition function is 

z = fv4> slow v4> fast eM-s slow -s% st -sl low -s% st ). (21) 

Thus I can integrate out Sf ast (\q\ > A/6 component) simply and obtain 

Z = Jvi> slow exp(-S slow -S 9 slmu ). (22) 

The remaining procedure of the renormalization is the scale transformation 

q — * q/b, w — * w/b and (j) ~^ 4>b 2 , (23) 

where I choose the dynamical exponent 1. The results are 

CO qO 

^slow ^ 

SLw - ^E E q 2 V(q/bMq,w)f 

w q=-A 

-A 

- s^E E q 2 v(q)\Hq,w)\\ (24) 

w q=—A 
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where I use V(q) — A ~ q 13 1 from the behaviours ([7]). Hence I obtain the renormalization group 
eq. 

4f = (X-fflf. (25) 
Substituting / = In 6 into this, I obtain renormalization group eqs. 

at f iv 

The TL parameter K is not renormalized but it shifts due to the constant A. 
P = l 

The dispersion relation of the Coulomb interaction includes the marginal part w ~ q and w ~ 



qy/\ lng| as well as > 1 case. Integrating out the fast moving part, I obtain the effective action 
of the slow part 



A/6 A/6 

Sslow = E E f (vq 2 +w 2 /vMq,w)\ 2 +gY, £ ?M<7M^)P 
w q=-A/b w q =-A/b 



(27) 



where I dare to leave the velocity in the Gaussian part. Note that I need not the renormalization 
of the velocity in the case /3 > 1. After the scale transformation, I obtain the eqs. 

I = 

oft 

d . v . 

5<^> - °> < 28 > 

where A is the constant appearing in ([7]). I see that K and the velocity v is renormalized instead 
of the no renormalization of g. The forward scattering become relevant through K, v and drive 
the system away from the TL fixed point. Note that this result holds irrespective of any filling 
kp. From these eqs., the size dependences of v and K are given by 

v(b) ~ Vl nL 



K(b) ~ 1/VlnL. (29) 
The velocity diverges weakly for long distances, which is consistent with the estimations of 



v 



^ from the behaviours d7J). 
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(3 = 3 

I use V(q) = A + Bq 2 lnq + Cq 2 H in the behaviours ([3). The g term of ([30]) is 

A/ft 

£ £ <Z 2 M0)g 2 lug + 5 2(0)g 2 ), (30) 

V) q=-\/b 

where the couplings <?i(0) and 52(0) are defined by gB and gC respectively. For the scale 
transformation (f23l) . the g term is changed to 

A 

£ £ q 2 [ 9 i(0)q 2 lnq/b 2 + (g 2 (0)/b 2 + gi (0)q 2 In q/b 2 ))]. (31) 

W q=~X 

Thus I obtain 

9i(b) = 1^(0) 

52(6) = 5i(0)plni + 52 (0)^. (32) 

-251 (0 

-2 52 (0-9i(0- (33) 



By / = /n6, I write this as 



dgijl) 
dl 

dg 2 (l) 
dl 

f3>3 

This case is same as eq. ([25]) putting f3 = 3. 



2. CFT in the TL liquid with LRI 

The Hamiltonian in the finite strip from the action ([3]) is 

H = H TL + g da 1 da 2 d ai (p(ai)d a2 (p(a2)V(\a 1 - a 2 \)0(\a 1 - a 2 \ - a ), (34) 
Jd 

where -Htl is TL liquid and D means the region D = {\a% — a 2 \ < L,—L/2 < a\,a 2 < L/2}. 
I introduce the step function 9(x) to avoid the ultra violet divergences which come from V(x) 
and the operator product expansion of d a <j)(a). For the small perturbation g the ground state 
energy E g varies as 

E g -Eg = g daid<T 2 V(\ai- a 2 \) <0\d (T1 (j){ai)d (T2 (j)(<T2)\0> 9{\ai- a 2 \ -ao) 
Jd 

daida 2 V(\ai - <r 2 |)[< Q\d Wl (p(wi)d W2 (p(w 2 )\Q > 
+ < Old^i^wi^^u^lO >] Tl=T2=Q 6(\ai - a 2 \ - a ), (35) 
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where I introduce the coordinates w = r + ia ( —L/2 < a < L/2, — oo < r < oo) and |0 > 
is the ground state of i?TL- From the characters of the Gaussian part (TL liquid part) I can 
separate as 4>(a,r) = ip(w) + <f(w) and derive < Q\dyj 1 (p{w\)d W2 ^p{w2)\0 >= 0. The content of 
the brackets is modified as follows: 

[< 0\d wl ip(wi)d W2 <p(w 2 )\0 > + < 1 ) 6»^ 2 <^ ( w 2 ) 1 >] 

Tl =T2 =0 

4 ll L J Zl (l_a)2 J ^(1-1)2^=^=0 



■K"/2tt 2 1 

~T ( T } 2 sin 2 ^r 2) ' (36) 

where I transform the correlation function < d Zl (p{zi)d Z2 (p{z2) >= ^ z ^_ Z2 yi in oo x oo z plane 
to that in the strip w thorough z = exp ^H. At present case d w (p{w) {du,(p(w)) have the spin 
s = 1(— 1) and conformal dimension A = 1(A = 1). Hence I obtain 

E'-E, = i^lqff 1 ' 2 ^ 1 (37) 

4 i/ 7-1/2 (sm7r|x |J" sm 7tx i/ 

where I impose the periodic boundary condition and use the interaction potential V(x) = 
l/(£ sin(^)) /? . Putting e = q /L for convenience, I give the differential of the integral part: 

d f 1 / 2 1 2 

7k / dx V | ■ 2 ^ (l^'l - £ ) = ~ 7 ■ I ■ 2 • (38) 

ue J-i/2 (sm7r|x l)^ sm 7rcc (sm7r|e|) p sin 7re 

After integrating the Taylor expansion about e of this quantity, I obtain 



f 1/2 , ' 1 1 „„ ', 2 r (vre)-^- 1 /3 + 2 . . fl+1 

+oa^r l,+ % (39) 

where (5 ^ odd. Therefore I can write the corrections in the form: 

4-* - f !^ + B(fflL + ^ + ^> +0 (^)]. (40) 

Here £?(/?), C{(3) and are given by 

-8-1 



1 + P 



c <« - w^" 1 (41> 

I can obtain A{0) by evaluating the above integral numerically. The result is shown in Fig. H^a). 
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For 0= odd, there exists the logarithmic correction instead of the eq. (|30]l . The results for 



respective are 



E-E n 



2[A +BL+ c ln i +D ^ + 0{ ^ )} for 13 = 1 

g[& + BL+C +jD ^lni + 0{^)\ for = 3 

<7[-^ + BL + £ + Z^+£-^lnI+0(-^)] for /? = 5 



(42) 



The C terms in eqs. (140p and (j42J) contribute to deviation of the central charge. The present LRI 
inevitably contains the contribution from the short range interaction: 5(x). The C term does 
not come from the short range types of interactions because the vacuum expected value ((d x 4>) 2 } 
vanishes. Thus C term is intrinsic in the present system with the LRI under periodic boundary 
condition. Because the velocity is not renormalized as I have seen from the renormalization 
group eqs., the C term contributes to the deviations of the effective central charge from the 
ground state energy. 

Next I derive the corrections for the energy of the excited state: 

E„ — E n 



g / daida 2 V(\ai - cr 2 |) < n\d ai <l)(a 1 )d a2 (j)(a2)\n > 9{\a\ - a 2 \ - a ) 
Jd 
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/ daida 2 V(\ai — a 2 \) } y < n\d ai (j)(ai)\a >< a|9 CT2 ^(cr 2 )|n > 9{\cr\ - cr 2 | 
Jd „ 



LP 



1/2 



dy 



(sin7r|y|)^ 



cos27r(s„ - s a )y 0(\y\ 



OLQ, 

L 



(43) 



where I use the results by Cardy [22J: 



< n\(b(a)\a > 



2ti\ 



C ■ {—\ x ip 



2-7ri(s n — s tt )g 



(44) 



Here j means dtp and \n > is the excited state of i?TL- 1 can derive the size dependence of 
eq. (|43p from likewise treatments as the ground state. After taking the derivative about 1/L, I 
expand about 1/L. Integrating them, I obtain 



E„ — Er, 



16TT 2 g C n j a C a j n { - " L/3 



,3) 



, BUS) , C{s n -s a ,f3) 

+ t - H t? — 



+ D(s n -s a ,p) +Q( ^ )} 



+ odd 



16tt 2 5 CnjaCajni ^^ + #L ln X + C ( S n ~ + D ( S n ~ S a )^U + O(^)} = 1 

16vr 2 5 J2 a C„ ja C ajn { Aisn L 7 a) +B{ + C{s n -s a )^\n{+ D(s n - s a )± + O(^)} = 3 
16vr 2 5 E Q C nja C ajn { A{Sn L - Sa) +B± + C{s n - s a )^ + D(s n - s a ) & In \ + 0{^)} = 5 



(45) 
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where B are the constant independent of s n ,s a . Here for /3 ^ odd, B{(5), C{s n — s a ,/3) and 
D(s n - s a ,f3) are given by 

B((3) ' 



(a 7r)^n(P - 1) 



C(s n -s a ,f3) = 1 [f - 2(s w - s a ) 2 } 1 (46) 



(aovr)^ 3 " 5 3 180 ; 72 3 J vr(/5-5)' 

It is not straightforward to determine A(s n — s a ,f3) generally. However about one particle 
excitation(s n = 0), I can obtain A(/3,s a ), which is shown for some s a in Fig. 2Kb) . Actually 
further consideration about the operator product expansion leads s n = s a = (see section 4.). 

I refer to the 0(1/L) dependences. These are due to the fact that the LRI includes the short 
range type interaction. Actually I can derive the same form 

~Y ^ ] C n j a C a j n (47) 

a 

as an ordinary finite scaling by replacing as V(|x|)0(|a;| — ao) — > S(x). As (deft) 2 is a part of the 
TL liquid, the 0(1/L) term can be erased under subtracting such the contributions first. Thus 
the 0(1/ L) term is not intrinsic. 

Summarizing the discussions in this appendix, I can prove that the Hamiltonian (|34p is 
described by c = 1 CFT for (3 > 1 in the excitation energy. However the effective central charge 
from the ground state depends on the interaction and deviates from 1. I find the nontrivial 
behaviors when (3 = odd, which corresponds to the integer points of the modified Bessel function 
as appearing in the behaviours ([7]). 
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Abstract 

I attempt to construct U(l) conformal field theory(CFT) in the Tomonaga-Luttinger(TL) 
liquid with 1 jr^ long-range interaction (LRI) . Treating the long-range forward scattering as 
a perturbation and applying CFT to it, I derive the finite size scalings which depend on the 
power of the LRI. The obtained finite size scalings give the nontrivial behaviours when (3 is 
odd and is close to 2. I find the consistency between the analytical arguments and numerical 
results in the finite size scaling of energy. 

1 Introduction 

Electron systems have attracted our much attention in the low energy physics. As the dimension 
of the electron systems decrease, the charge screening effects become less important. In spite of 
these facts, models with short-range interaction have been adopted in many researches of one 
dimensional electron systems. The recent advanced technology makes it possible to fabricate 
quasi-one-dimension systems. Actually in low temperature the effect of Coulomb force has been 
observed in GaAs quantum wires [Tj , quasi-one-dimensional conductors O El H] and ID Carbon 
nanotubes 0I51E]. 

The systems with 1/r Coulomb repulsive forward scattering was investigated on the long 
distance properties by bosonization techniques [SJ. The charge correlation function decays with 
the distance as exp(— const. (In x) 1 / 2 ) more slowly than any power law. The momentum distri- 
bution function and the density of state does not show the simple power law singular behaviour. 
The logarithmic behaviours appear in the power [9]. These mean that the system is driven to 
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the Wigner crystal which is quite different from the ordinary TL liquid. The investigation for 
the interaction l/r 1 ~ e through the path integral approach [10] reconfirms the slower decaying of 
the single particle Green function for e = and leads the faster decay for < e(<C 1) than any 
power type. 

The numerical calculation in the electron system with the Coulomb interaction shows that 
the larger range of the interaction causes the insulator (charge density wave) to metal (metallic 
Wigner crystal) transition [11]. In the spinless fermion system, the convergence of the Luttinger 
parameters exhibits the quasi-metallic behaviour different from the simple TL one |12j . 

As I will discuss below, the forward scattering is irrelevant for (3 > 1. As an instance of the 
effect of the long-range Umklapp scattering, it was reported that the 1/r 2 interaction makes the 
system gapless to gapful through the generalized Kosterlitz-Thouless transition [13j . 

In this paper I discuss CFT in the system with LRI. The basic assumptions of CFT are 
symmetries of translation, rotation, scale and special conformal transformation. Besides them I 
assume the short-range interaction in the CFT. Hence it is a subtle problem whether the CFT 
can describe the system with LRI. 

Of LRIs, up to now, the solvable models with 1/r 2 interaction were discussed |14 4 115 1 [TB I 117]. 
With the Bethe ansatz, the conformal anomaly and the conformal dimensions were calculated 
and the system proved to be described by c = 1 CFT. In fact the central charge from the specific 
heat agrees with c = 1. On the other hand, the ground state energy is affected by the LRI and 
the periodic nature. The effective central charged deviates from c = 1. 

In general, the CFT for LRI which breaks the locality, has been left as unsettled problem. It 
is significant to clarify the validity of the CFT to the systems with LRI. I investigate the tight- 
binding model with l/r@ interaction as one of such problems. The low energy effective model 
consists of TL liquid, the long-range forward scattering and the long-range spatially oscillating 
Umklapp scattering. Extending arguments appearing in Ref. [18] to the TL liquid with the 
long-range forward scattering, I derive the finite size scalings. In the tight-binding model with 
\j T P interaction, I calculate numerically the size dependences of energy and the coefficients of 
1/L y . And I see numerically the relations between the velocity, susceptibility and Drude weight, 
which CFT requires. 

2 Field theoretical approach 

I consider the following tight-binding Hamiltonian of the interacting spinless Fermions: 

L L 

h = - j>W + m + 1 2> - uma - *> - 1/2), w 

X I define the "effective central charge" by c = — in the finite size scaling of the ground state energy E g — 
aL — . I use the word "effective central charge" in this sense. 
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where the operator cj (cJ) annihilates (creates) the spinless Fermion in the site j and pj = CjCj 
is the density operator. In order to treat this model under the periodic boundary condition, 
I define the chord distance between the sites i and j: r%j = (— sin nfcj) ) where L is the site 
number. Using this, I express the LRI as V(i — j) = — T ^u-iW „ • 

(_ sin " L JU )P 

By the bosonization technique, I obtain the effective action of the Hamiltonian (pQ) for the 
arbitrary filling: 

S = drdx — (V(j)) 2 +g / drdxdx d x (f>(x , r)V (x — x )di<j)(x ,r) 

J 27TA J 

+ g j drdxdx cos(2fcpx + V2<f>(x, r))V(x — x ) cos(2fcpx + V2<ft(x , t)) , (2) 

where V(x) = K is the TL parameter and kF is the Fermi wave number. And g is the 
couping constant proportional to g. The first term of ([2]) is the TL liquid and the second term 
is the long-range forward scattering. The last term is the spatially oscillating Umklapp process 
which includes cos 2y2(f) which comes from the interaction between the neighbour sites. 

Schulz analyzed the effects of the Coulomb forward scattering by the bosonization technique 
in the electron system|8j. He discussed the quasi- Wigner crystal of electrons due to the Coulomb 
forward scattering. Here I focus on the effects of the 1/r 13 forward scattering in the spinless 
Fermions system. I treat the action : 

drdx- — — (V0) 2 + g j drdxdx d x (f>(x, t)V(\x — x \)d x /<p(x , r) (3) 



2-kK 

for any filling hp. To investigate in the Fourier space, I choose the form V(x) = ^ 2+ ^ 2 ^/ 2 , 
where a is the ultra-violet cut-off. In the wave number space, the action ([3]) is expressed as 

/2ir 
dqdw{-(q 2 + w 2 ) + gq 2 V Xq)}\<Kq, w)\ 2 , (4) 

where V{q) is the Fourier transformation of V(x): 

n«) = r( ^ /2 2 ) ^ F / 2 -i /2 (^) /3/2 ' 1/2 ^-v^ a g)- ( 5 ) 

Here K u {x) is the modified Bessel function of z^th order and T(x) is the gamma function. From 
this, the dispersion relation is 

W " = 9 2 {1 + ^^)}- (6) 
The long wavelength behaviors of V(q) are given by 

A + B(qa) 2 + C{qaf- 1 + ■■■ (3 > and f3 + odd 

A + B In ga H = 1 

V(q) ~ { A + B{qa) 2 \nqa + C{qa) 2 + ■■■ (3 = 3 (7) 

A + B(qa) 2 + C{qaf ln(ga) + D{qa) A + ■■■ (3 = 5 
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where A, B, C and D are the functions of j3. For the case where (3 > and ^ odd, the coefficient 
B = B((3), C = C{p) is given by 

vr 3 / 2 f 
B{J3 ^ ~ 4 2 / 3 / 2 - 1 /2r(^)r(/3/2)sin^^ 

c <« - -' 3/2 ^- lr( m )r ; /2)sin i^ ' < 8 > 

From the eqs. ([6]) and ([7]), I see that (</a)^ _1 and In go terms for < (3 < 1 affect the linear 
dispersion essentially. Especially for (3 = 1, there is the analysis by Schulz, where the charge 
density correlation function is calculated [8]. According to it, in the present spinless case, the 
LRI drives the ground state from the TL liquid to the quasi- Wigner crystal as (3 — > 1+. The 
slowest decaying part of the density correlation function is given by 



< p(x)p(0) > ~ cos(2/cpx)exp(— cy^og x), (9) 

where c is a function of K, which exhibits slower spatial decay than the power decay of TL 
liquid. 

Then I see the effects of the long-range forward scattering in the standpoints of the renor- 
malization of g. The renormalization group eqs. of g, v and K are simply derived for the 
long wave-length (see Appendixes.). From the renormalization eqs., the g terms are relevant for 
(3 < 1, marginal for (3 = 1 and irrelevant for (3 > 1. Thus it is expected that the system becomes 
the quasi- Wigner crystal caused by the forward scattering for (3 < 1 and the system becomes 
the TL liquid when (3 > 1. I see that 

= J dx'd x <j)(x,T)V(x -x')d x ><t>(x ,t) (10) 

has the scaling dimension x g = (3 + 1 for 1 < (3 < 3 and 4 for (3 > 3. As the weak logarithmic 
corrections appear for /3=odd, I here distinguish <3?(x) for /3=3 from the scaling functions. I also 
find the consistency on these scaling dimensions by CFT. By using the first order perturbation, I 
can know the effects of the long-range forward scattering. Based on CFT, the finite size scalings 
of energies for no perturbations are given [191 1201 EH [22] by 

„ ^ 2irvx n 

= — 

Eg = e 9 L -jr> (11) 
where x n is the scaling dimension of the primary field denoted by n, v is the sound velocity and 
c is the central charge. Considering the LRI, I can extract the corrections to these energy size 
scalings(see Appendixes.): 

AE n = ^(l + ^^ + 0(l/L 2 )) 
E g = ( eg+9 (0)const.)L-— (c + 5 (0) const. +5(0)-^^ + 0(1/L 2 )), 

(12) 
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where /?(> 1) is not odd. And the constants are the functions of (3. Note that for (3 = odd 
cases, the logarithmic corrections appear. They correspond to the integer points of the modified 
Bessel function, which appear in the long- wave behaviours ©• I can reproduce these anomalies 
for (3 = odd by the CFT. Moreover from CFT I can show that there are the anomalies in the 
general excitations and the ground state energy. The details are shown in Appendixes. The 
0(1/L 2 ) terms come from the irrelevant field L_2-^-2l and the long-range g term. The first eq. 
of (|12p means that the long-range forward scattering has the scaling dimension x g = (3 + 1 
for 1 < (3 < 3 and 4 for (3 > 3 effectively. These respective scaling dimensions are consistent 
with the estimation from the renormalization group eqs. of g, that I mentioned above (see 
Appendixes.). 

The energy finite size scalings (|12p mean that the LRI has the higher order influences than 
1/L to the excitation energy and the LRI affects the 1/L term in the finite size scaling of the 
ground state energy. Here I note that it becomes difficult to calculate the central charge from 
finite size scalings ([lip unless the effects of the LRI to 0(1/ L) terms are known. 

It is notable to compare the eqs. ()12p with the case where the perturbations are of short- 
range type. Ludwig and Cardy calculated the contributions of the short-range perturbation |18j . 
The results for the irrelevant perturbation, — g^2 r 4>(r), which has the scaling dimensions x > 2 
are 



A En = ^(i + ^cw^r 2 ) 

j-j Xfi j-j 

7T11 

E g = (e g + g^) CO nst.)L-—(c + g(Qf^= + 0(l/L^)), (13) 



7™/ /^2 COI1S t- /rlT-fi\ 



where the O(g) terms do not appear in the ground state scaling because we set (<p) = for 
the short-range interaction. These results mean that the x > 2 irrelevant field has influences of 
the higher order to the finite size scalings (jlip . And their result contains parts not so simple. 
There are the special points of scaling dimension x = 1,3,5, and x = 2 which is related to the 
appearance of logarithmic corrections. 

To the contrary, I see ($) ^ in the long-range case, where $ is defined in eq. (jlOp . The 
LRI gives the 0(\/L) intrinsic influence to the finite size scaling of the ground state energy, as 
appearing in the scalings (|12p . even if the LRI is irrelevant, that is, x g > 2. 

3 Numerical calculations 

Through the Jordan- Wigner transformation, I transform the model ([1]) to S = 1/2 spin Hamil- 
tonian for the numerical calculations: 

H = -E( s / 5 7+i +h - c ) + fE 5 ^-^i) 5 i- ( 14 ) 

3 i+3 
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I impose the periodic boundary condition S^+i = Si to this model. Using the Lanczos algorithm 
I perform the numerical calculations for the Hamiltonian (|14p . 

I have found analytically the corrections to the energy scalings (|lip caused by the long- 
range forward scattering. If the oscillating Umklapp process term of ([2]) is irrelevant and does 
not disturb the energy scalings, the finite size corrections due to the forward scattering are 
expected to appear in the excited state energies and the ground state energy. I attempt to 
detect the contribution of the forward scattering. 

I numerically calculate the size dependences of the excitation energy AE(m = 1/L) and the 
ground state energy E g (m = 0), E g (m = 1/L) for g = 0.5. Here I define the magnetization 
m = Y2j Sj/L w hich is the conserved quantity. Fitting the one particle excitation energy as 
LAE(m = 1/L) = a+-£s + jjz,l show the power c versus the powers j3 in Fig. [TJ I see the power 
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Figure 1: The numerically calculated powers c in the excitation energy LAE(m = 1/L) are 
shown versus (3 for g = 0.5. Here I use the scaling form: LAE = a + + where a, b, c and 
d are determined numerically. If the LRI is not present, the energy finite size scaling must take 
the form :LAE = A + -^ + 0(1/L 4 ), where A and B are constant values. 



c agrees with theoretical predictions: (3 — 1 except for (5 = 2. 1 shall discuss the (3 = 2 case later. 
Fitting the ground state energy per site as E g /L = a + jp + , I plot the powers d versus (3 
in Fig. [2j I see that the power d do not show agreements with theoretical predictions (3 + 1 in 
E g (m = 0)/L. These disagreements may be caused by the oscillating Umklapp process which 
becomes relevant at only m = filling. On the contrary, the oscillating Umklapp process is 
irrelevant at m / 0. Actually, in Fig. [51 I see that the power d show agreements with theoretical 
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Figure 2: The numerically calculated powers d in the ground state energies E g {m = 1/L)/L and 
E g (m = 0) / L are shown versus (3 for 51 = 0.5. Here I use the scaling form: E g /L = a + + 
where o, 6, c and d are determined numerically. If the LRI is not present, the energy finite size 
scaling must take the form :E g /L = A + + j^, where A, B and C are constant values. 
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predictions (3+1 in E g (m = 1/L)/L except for (5 = 2. 

As I stated above, for (3 = 2, the power c in the excitation energy LAE(m = 1/L) = 
a + jc + -rz apparently shows disagreement with theoretical value (3—1 and likewise for (3 = 2, 
the power d in the ground state energy E g (m = 1/L)/L = a + + jd apparently shows 
disagreement with theoretical value (3 + 1. I investigate the reason for these disagreements. 
In Fig. [3] I show the numerically obtained coefficient of 1/L d in the size scalings E g (ra = 
0)/L, E g {m = 1/L)/L = a+ + and the numerically obtained coefficient of 1/L C in the size 
scaling LAE(m = l/L)=a + j^ + jj I . I observe that the coefficient of 1/L C in LAE(m = 1/L) 
and the coefficient of 1/L d in E g (m = 1/L)/L become small around (3 = 2. So for (3 = 2, 
1/L 2 dependence appears rather than 1/L in LAEim = 1/L) (see Fig. [TJ). Likely for (3 = 2, 
1/L 4 dependence appears rather than 1/L 3 in E g (m = 1/L)/L (see Fig. [2j). I observe that the 
coefficient of 1/L d in E g (m = 0)/L show the different behaviour from those in E g (m = 1/L)/L 
in Fig. This difference may come from the spatially oscillating Umklapp process that opens 
the gap at m = and disturbs the finite size scaling. 



H — ' 

£= 
Q 

O 
O 



10 

8 
6 
4 
2 

-2 
-4 



in E g (m=0)/L 
in E g (m=l/L)/L 
in LAE(m= 1/L) 





1.5 



2.5 



|6 



3.5 



Figure 3: I show the numerically obtained coefficients of 1/L d in the size scalings E g (m = 
0)/L, E g (m = 1/L)/L = a + + j-& and the numerically obtained coefficient of 1/L C in the size 
scaling LAE(m = l/L) = a+^ + -^. I observe that the coefficient of 1/L d in E g (m = 1/L)/L 
and the coefficient of 1/L C in LAE(m = 1/L) become small around (3 = 2. The coefficients of 
1/L d in E g {m = 0)/L show the different behaviour from that in E g {m = 1/ L)/L. This difference 
may be caused by the spatially oscillating Umklapp process term. 



I can obtain A{(3) in the scalings (|40p and (|45p by evaluating the integrals. The results are 
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shown in Fig. 0] (a) and (b). The analytical A{(5) in the scalings (|40p and A(f3, s) in the scalings 
(|45p for s = fit with the points in Fig. [3J well. The curve only for s = in Fig. IHb) shows 
the good fitting. This point shall be discussed later. These reveal that the present numerical 
calculation of the tight-binding model agrees with the CFT analysis of the long-range forward 
scattering. 




3.5 




Figure 4: (a) A(P), the coefficient of in the eq. (I40p is shown. I see that A(j3) has zero 

point close to j3 = 2. This curve coincides with the results from the numerical calculation in the 
tight-binding model shown in Fig. [3j (b) A((3), the coefficient of 1/L@, in the eq. (|45p is shown 
for some s. Analytically only s = is meaningful for particle excitations. A(f3) for s = has 
zero point close to (3 = 2. This coincides with the results from the numerical calculation in the 
tight-binding model shown in Fig. [31 

Next I survey whether the long-range tight-binding model satisfies the necessary condition of 
CFT. The operator cos V2<ft has the scaling dimensions K/2 and the operator e =tl v / 2 e h as 
in the regime of the TL liquid. The two quantities 2K/v and vK/2 are the compressibility and 
the Drude weight respectively in the regime of the TL liquid. If c = 1 CFT is valid to the 
tight-binding model with the LRI, the two quantities are related to the two excitations with the 
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symmetries q = ir, m = and q = ir,m = 1/L respectively: 

2K/v = l/(LAE(m=l/L,q = ir))=x 

vK/2 = LAE(q = ir)=D. (15) 

I show the numerically calculated quantities \ and D in Fig. [5] and [U where I use the sizes L = 
16, 18 and 20 and extrapolate the data. For g < 0, x (which is the susceptibility, irrespective 
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Figure 5: The extrapolated K/v(= x/2) is plotted versus the strength g. I use the scaling form 
v/K{L) = v/K(oo) + jjhr for < 3, and v/K{L) = v/K{oo) + ^ for p > 3, where v/K(oo), 
a is determined numerically. 

of the CFT arguments) exhibits the rapid increase which suggests the phase separation. In spin 
variables' language for (pQ), this phase separation is nothing but the ferromagnetic phase. Hence 
for the larger (3 the point of the phase separation approaches to —1. For g > I see the weak 
tendency that the the quantity x becomes smaller as (3 is smaller for g less than about 1. I find 
that the the quantity D of g > become larger as (3 approaches to (3 = 1. 

In Fig. [7]I plot the velocity versus the strength g for the various powers (3, where the velocity 
is defined by 



^-AE(q = 27T/L). 



(16) 



I see that the velocities are finite values for f3 > 1 , as is expected. There are the points where 
the velocities are zero, implying the phase separation. 

In Fig. [8] I plot the quantity versus the strength g for the various powers f3. If the 
present system is described by c = 1 CFT, this quantity is 1 from eqs. (|15p . I find the regions 
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Figure 6: The extrapolated vK(= 2D) is plotted versus the strength g. I use the scaling form 
LAE = a + -jj, where a,b and c are determined numerically. 




g 



Figure 7: The extrapolated spin wave velocity v is plotted versus the strength g. I use the 
scaling form LAE = a + where a,b and c are determined numerically. 
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where = 1 in Fig. [gj The regions become wider as (3 approaches to 1 for g > 0. For larger 
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Figure 8: The normalization is plotted versus the strength g. 



g, the normalization breaks owing to the generations of mass. 



4 Discussion 



I have investigated the system with the 1/r 13 interaction by applying CFT to it and by the 
numerical calculation. At first I have analyzed TL liquid with the l/r@ forward scattering by 
utilizing the CFT and I have found that the l/r@ forward scattering works as higher order 
corrections in the excitation energy, whereas the effective central charge in the scaling of the 
ground state energy depends on the interaction and it deviates from c = 1. The deviation are 
like the solvable 1/r 2 models |14| [15} [TBI fT7] . Next I have numerically calculated the ground 
state energy and excitations energies in the tight-binding model with l/r@ interaction, which is 
expected to include the above \/r@ forward scattering in the low energy. The numerical results 
are in accordance with the analysis with CFT of the long-range forward scattering. Furthermore 
I have numerically checked the normalization -^r = 1, which is the necessary condition for c = 1 



CFT. 



For (3 ~ 2, the coefficient A{f3) in the ground state energy vanishes. This seems to correspond 
to the exact solution for (3 = 2 [17] which states that the finite size scaling of ground state has 
no higher order term than 1/L. The coefficient D((3) of 1/L 3 in eq. (I42p does not vanish for 
(3 = 2. However the present argument is the first order perturbation theory. With higher order 
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treatments, I may clarify this. In any case, with consistency in many points I could construct 
CFT in the system with non-local interaction. 

The numerical calculations in the tight-binding model support the finite size scalings (|40|) 
and (|45p . In one particle excitation energy LAE(m = l/L), the coefficients of l/L^ fit with 
s = case in Fig. SJ The coefficients from the long-range forward scattering are related with 
the operator product expansion. I can prove that only s = case is relevant for the particle 
excitation. Using < ip(z)(p(z') >= — ln(z — z) and < (p(z)(p{z) >= — -f- ln(z — z'), I confirm 
the operator product expansions: 

dtp{z):<*&^: = 1 , : e ^(*W) . +reg . 

4 z — z 

3<p(z) : ^ '>*') : = ^JL. : e*^<-V) . +reg . (17) 



4 


z — z 


iv^ 


1 


4 f 


=- '. 

1 

— z 


1 


1 


4K (z 


- z'y 


1 


1 


4K (z 


- z'y 



T(z) : e *^(,V) : = : e^V) : : e ^*(*V) :) : +reg. 

iv z — z 

T(z) : e ^(-'- ; ') : = -L 1 : e ^(,',;') . . e ^(,', 2 "') :) . +reg _ ; 

1 '' ' ' K z — z 

where I define T(z) = -j?(d<p(z)) 2 , f{z) = --|(^(z)) 2 and = j?((p(z) - <p{z)). From 

the first and the second eqs., I see C a \o = —iV2/4, C al Q = iy/2/4 for a = 1 and C a \Q = C al Q = 
otherwise, where 0( ) and 1 denote d<p(z) ( d(p{z)) ) and : e l ^ s ( z ' z ) : . From the third and 
the fourth eqs., I see : e l ^ e ( z ^ z ) • have the conformal dimension (1/4K, 1/4K) and spin 0. As 
i(dip(z) — d(p{z))/2 is associated with d a {a) for z = exp(^^), I obtain 



< a\d a (a)\l > 



^i(C 110 -C m )/2 = ^ fora = l 
otherwise, 



which means that only s = is relevant for the particle excitation and the last eq. in (|43p has 
no cosine term. 

I would discuss the size effects for (3 = 1. As seen in the eqs. (|29|) and (|45|) . the velocity shows 
the weak divergence for the size and the Luttinger parameter vanishes gradually for increasing 
size. This is consistent with the numerical tendency (see Figs. 8 and 9 in Ref. [E]). The size 
effect of the Drude weight (proportional to the charge stiffness) is now given by vk/2 ~ const, 
as the logarithmic contributions cancel. The numerical data (see Fig. 7 in Ref. [12]) shows 
the metallic behaviour at small and intermediate magnitude interaction strength (larger than 
the CDW transition point V = 2 by the short range interaction). I think that the long-range 
forward scattering enhances the metallic character. For fairy large interaction the long-range 
Umklapp scattering becomes relevant and the charge stiffness is suppressed. Finally I would 
add the size effect of the compressibility: 

LAE(p = l/2 + l/L) = 1/x = OQnL) -> oo, (19) 
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which comes from the results (|29[) by the RG analysis and the CFT arguments. The compress- 
ibility x goes to weakly for increasing size. 

To summarize, within the perturbation theory I have constructed CFT in the TL liquid 
with l/r@ long-range forward scattering. I have found that the interaction gives the nontrivial 
behaviour for (3 =odd and j3 ~ 2. I have numerically checked the finite size scalings obtained 
from CFT in the tight-binding model with l/r@ LRI. Our analysis and numerical calculations 
exhibit consistency with each other. 

Appendixes 

1. Renormalization group equation 

0</3<lorl</3<3 

I derive the renormalization group equations heuristically. Let us start from the action (J3|): 

S = E E §(9 2 +™ 2 M^)| 2 +5£ E q 2 V(q)\<t>(q,w)\ 2 
w q=—A w q=—A 

A/6 -A/6 A A/6 -A/6 A 

= Ei E + E + Ei+^E{ E + E + E)- ( 2 °) 

w q=-A/b <j=-A g=A/6 w q =-A/b q=-A q=A/b 

The partition function is 

z = fv4> slow v4> fast eM-s slow -s% st -sl low -s% st ). (21) 

Thus I can integrate out Sf ast (\q\ > A/6 component) simply and obtain 

Z = Jvi> slow exp(-S slow -S 9 slmu ). (22) 

The remaining procedure of the renormalization is the scale transformation 

q — * q/b, w — * w/b and (j) ~^ 4>b 2 , (23) 

where I choose the dynamical exponent 1. The results are 

CO qO 

^slow ^ 

SLw - ^E E q 2 V(q/bMq,w)f 

w q=-A 

-A 

- s^E E q 2 v(q)\Hq,w)\\ (24) 

w q=—A 
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where I use V(q) — A ~ q 13 1 from the behaviours ([7]). Hence I obtain the renormalization group 
eq. 

4f = (X-fflf. (25) 
Substituting / = In 6 into this, I obtain renormalization group eqs. 

at f iv 

The TL parameter K is not renormalized but it shifts due to the constant A. 
P = l 

The dispersion relation of the Coulomb interaction includes the marginal part w ~ q and w ~ 



qy/\ lng| as well as > 1 case. Integrating out the fast moving part, I obtain the effective action 
of the slow part 



A/6 A/6 

Sslow = E E f (vq 2 +w 2 /vMq,w)\ 2 +gY, £ ?M<7M^)P 
w q=-A/b w q =-A/b 



(27) 



where I dare to leave the velocity in the Gaussian part. Note that I need not the renormalization 
of the velocity in the case /3 > 1. After the scale transformation, I obtain the eqs. 

I = 

oft 

d . v . 

5<^> - °> < 28 > 

where A is the constant appearing in ([7]). I see that K and the velocity v is renormalized instead 
of the no renormalization of g. The forward scattering become relevant through K, v and drive 
the system away from the TL fixed point. Note that this result holds irrespective of any filling 
kp. From these eqs., the size dependences of v and K are given by 

v(b) ~ Vl nL 



K(b) ~ 1/VlnL. (29) 
The velocity diverges weakly for long distances, which is consistent with the estimations of 



v 



^ from the behaviours d7J). 
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(3 = 3 

I use V(q) = A + Bq 2 lnq + Cq 2 H in the behaviours ([3). The g term of ([30]) is 

A/ft 

£ £ <Z 2 M0)g 2 lug + 5 2(0)g 2 ), (30) 

V) q=-\/b 

where the couplings <?i(0) and 52(0) are defined by gB and gC respectively. For the scale 
transformation (f23l) . the g term is changed to 

A 

£ £ q 2 [ 9 i(0)q 2 lnq/b 2 + (g 2 (0)/b 2 + gi (0)q 2 In q/b 2 ))]. (31) 

W q=~X 

Thus I obtain 

9i(b) = 1^(0) 

52(6) = 5i(0)plni + 52 (0)^. (32) 

-251 (0 

-2 52 (0-9i(0- (33) 



By / = /n6, I write this as 



dgijl) 
dl 

dg 2 (l) 
dl 

f3>3 

This case is same as eq. ([25]) putting f3 = 3. 



2. CFT in the TL liquid with LRI 

The Hamiltonian in the finite strip from the action ([3]) is 

H = H TL + g da 1 da 2 d ai (p(ai)d a2 (p(a2)V(\a 1 - a 2 \)0(\a 1 - a 2 \ - a ), (34) 
Jd 

where -Htl is TL liquid and D means the region D = {\a% — a 2 \ < L,—L/2 < a\,a 2 < L/2}. 
I introduce the step function 9(x) to avoid the ultra violet divergences which come from V(x) 
and the operator product expansion of d a <j)(a). For the small perturbation g the ground state 
energy E g varies as 

E g -Eg = g daid<T 2 V(\ai- a 2 \) <0\d (T1 (j){ai)d (T2 (j)(<T2)\0> 9{\ai- a 2 \ -ao) 
Jd 

daida 2 V(\ai - <r 2 |)[< Q\d Wl (p(wi)d W2 (p(w 2 )\Q > 
+ < Old^i^wi^^u^lO >] Tl=T2=Q 6(\ai - a 2 \ - a ), (35) 
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where I introduce the coordinates w = r + ia ( —L/2 < a < L/2, — oo < r < oo) and |0 > 
is the ground state of i?TL- From the characters of the Gaussian part (TL liquid part) I can 
separate as 4>(a,r) = ip(w) + <f(w) and derive < Q\dyj 1 (p{w\)d W2 ^p{w2)\0 >= 0. The content of 
the brackets is modified as follows: 

[< 0\d wl ip(wi)d W2 <p(w 2 )\0 > + < 1 ) 6»^ 2 <^ ( w 2 ) 1 >] 

Tl =T2 =0 

4 ll L J Zl (l_a)2 J ^(1-1)2^=^=0 



■K"/2tt 2 1 

~T ( T } 2 sin 2 ^r 2) ' (36) 

where I transform the correlation function < d Zl (p{zi)d Z2 (p{z2) >= ^ z ^_ Z2 yi in oo x oo z plane 
to that in the strip w thorough z = exp ^H. At present case d w (p{w) {du,(p(w)) have the spin 
s = 1(— 1) and conformal dimension A = 1(A = 1). Hence I obtain 

E'-E, = i^lqff 1 ' 2 ^ 1 (37) 

4 i/ 7-1/2 (sm7r|x |J" sm 7tx i/ 

where I impose the periodic boundary condition and use the interaction potential V(x) = 
l/(£ sin(^)) /? . Putting e = q /L for convenience, I give the differential of the integral part: 

d f 1 / 2 1 2 

7k / dx V | ■ 2 ^ (l^'l - £ ) = ~ 7 ■ I ■ 2 • (38) 

ue J-i/2 (sm7r|x l)^ sm 7rcc (sm7r|e|) p sin 7re 

After integrating the Taylor expansion about e of this quantity, I obtain 



f 1/2 , ' 1 1 „„ ', 2 r (vre)-^- 1 /3 + 2 . . fl+1 

+oa^r l,+ % (39) 

where (5 ^ odd. Therefore I can write the corrections in the form: 

4-* - f !^ + B(fflL + ^ + ^> +0 (^)]. (40) 

Here £?(/?), C{(3) and are given by 

-8-1 



1 + P 



c <« - w^" 1 (41> 

I can obtain A{0) by evaluating the above integral numerically. The result is shown in Fig. H^a). 
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For 0= odd, there exists the logarithmic correction instead of the eq. (|30]l . The results for 



respective are 



E-E n 



2[A +BL+ c ln i +D ^ + 0{ ^ )} for 13 = 1 

g[& + BL+C +jD ^lni + 0{^)\ for = 3 

<7[-^ + BL + £ + Z^+£-^lnI+0(-^)] for /? = 5 



(42) 



The C terms in eqs. (140p and (j42J) contribute to deviation of the central charge. The present LRI 
inevitably contains the contribution from the short range interaction: 5(x). The C term does 
not come from the short range types of interactions because the vacuum expected value ((d x 4>) 2 } 
vanishes. Thus C term is intrinsic in the present system with the LRI under periodic boundary 
condition. Because the velocity is not renormalized as I have seen from the renormalization 
group eqs., the C term contributes to the deviations of the effective central charge from the 
ground state energy. 

Next I derive the corrections for the energy of the excited state: 

E„ — E n 



g / daida 2 V(\ai - cr 2 |) < n\d ai <l)(a 1 )d a2 (j)(a2)\n > 9{\a\ - a 2 \ - a ) 
Jd 
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/ daida 2 V(\ai — a 2 \) } y < n\d ai (j)(ai)\a >< a|9 CT2 ^(cr 2 )|n > 9{\cr\ - cr 2 | 
Jd „ 



LP 



1/2 



dy 



(sin7r|y|)^ 



cos27r(s„ - s a )y 0(\y\ 



OLQ, 

L 



(43) 



where I use the results by Cardy [22J: 



< n\(b(a)\a > 



2ti\ 



C ■ {—\ x ip 



2-7ri(s n — s tt )g 



(44) 



Here j means dtp and \n > is the excited state of i?TL- 1 can derive the size dependence of 
eq. (|43p from likewise treatments as the ground state. After taking the derivative about 1/L, I 
expand about 1/L. Integrating them, I obtain 



E„ — Er, 



16TT 2 g C n j a C a j n { - " L/3 



,3) 



, BUS) , C{s n -s a ,f3) 

+ t - H t? — 



+ D(s n -s a ,p) +Q( ^ )} 



+ odd 



16tt 2 5 CnjaCajni ^^ + #L ln X + C ( S n ~ + D ( S n ~ S a )^U + O(^)} = 1 

16vr 2 5 J2 a C„ ja C ajn { Aisn L 7 a) +B{ + C{s n -s a )^\n{+ D(s n - s a )± + O(^)} = 3 
16vr 2 5 E Q C nja C ajn { A{Sn L - Sa) +B± + C{s n - s a )^ + D(s n - s a ) & In \ + 0{^)} = 5 



(45) 



18 



where B are the constant independent of s n ,s a . Here for /3 ^ odd, B{(5), C{s n — s a ,/3) and 
D(s n - s a ,f3) are given by 

B((3) ' 



(a 7r)^n(P - 1) 



C(s n -s a ,f3) = 1 [f - 2(s w - s a ) 2 } 1 (46) 



(aovr)^ 3 " 5 3 180 ; 72 3 J vr(/5-5)' 

It is not straightforward to determine A(s n — s a ,f3) generally. However about one particle 
excitation(s n = 0), I can obtain A(/3,s a ), which is shown for some s a in Fig. 2Kb) . Actually 
further consideration about the operator product expansion leads s n = s a = (see section 4.). 

I refer to the 0(1/L) dependences. These are due to the fact that the LRI includes the short 
range type interaction. Actually I can derive the same form 

~Y ^ ] C n j a C a j n (47) 

a 

as an ordinary finite scaling by replacing as V(|x|)0(|a;| — ao) — > S(x). As (deft) 2 is a part of the 
TL liquid, the 0(1/L) term can be erased under subtracting such the contributions first. Thus 
the 0(1/ L) term is not intrinsic. 

Summarizing the discussions in this appendix, I can prove that the Hamiltonian (|34p is 
described by c = 1 CFT for (3 > 1 in the excitation energy. However the effective central charge 
from the ground state depends on the interaction and deviates from 1. I find the nontrivial 
behaviors when (3 = odd, which corresponds to the integer points of the modified Bessel function 
as appearing in the behaviours ([7]). 
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